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ABSTRACT 

We explore the detailed and broad properties of carbon burning in Super Asymptotic Giant Branch 
(SAGB) stars with 2755 MESA stellar evolution models. The location of first carbon ignition, quench¬ 
ing location of the carbon burning flames and flashes, angular frequency of the carbon core, and carbon 
core mass are studied as a function of the ZAMS mass, initial rotation rate, and mixing parameters 
such as convective overshoot, semiconvection, thermohaline and angular momentum transport. In gen¬ 
eral terms, we find these properties of carbon burning in SAGB models are not a strong function of the 
initial rotation profile, but are a sensitive function of the overshoot parameter. We quasi-analytically 
derive an approximate ignition density, pi gn ^ 2.1 x 10 6 g cm -3 , to predict the location of first car¬ 
bon ignition in models that ignite carbon off-center. We also find that overshoot moves the ZAMS 
mass boundaries where off-center carbon ignition occurs at a nearly uniform rate of AMzams/^/ov~ 
1.6 M 0 . For zero overshoot, / ov =0.0, our models in the ZAMS mass range ~ 8.9 to 11 M 0 show 
off-center carbon ignition. For canonical amounts of overshooting, / ov =0.016, the off-center carbon 
ignition range shifts to « 7.2 to 8.8 M 0 . Only systems with / ov > 0.01 and ZAMS mass ~ 7.2-8.0 
M 0 show carbon burning is quenched a significant distance from the center. These results suggest a 
careful assessment of overshoot modeling approximations on claims that carbon burning quenches an 
appreciable distance from the center of the carbon core. 

Subject headings: stars: evolution — stars: interiors — stars: rotation — supernovae: general 


1. INTRODUCTION 

When a single star on the main sequence (MS) ex¬ 
hausts the supply of hydrogen in its core, the core con¬ 
tracts and its temperature increases, while the outer lay¬ 
ers of the star expand and cool. The star becomes a red 
giant (e.g., Iben 1991; Stancliffe et al. 2009; Karakas & 
Lattanzio 2014). The subsequent onset of helium burn¬ 
ing in the core causes the star to populate the horizontal 
branch for more metal-poor stars or the red clump for 
more metal-rich stars (Cannon 1970; Faulkner & Can¬ 
non 1973; Seidel et al. 1987; Castellani et al. 1992; Girardi 
1999). After the star exhausts the supply of helium in its 
core, the carbon-oxygen (henceforth CO) core contracts 
while the envelope once again expands and cools along 
a path that is aligned with its previous red-giant track. 
The star becomes an asymptotic giant branch (AGB) 
star (e.g., Hansen et al. 2004; Herwig 2005; Kippenhahn 
et al. 2012; Salaris et al. 2014; Fishlock et al. 2014). 

The minimum mass for carbon ignition is usually re¬ 
ferred to as M up 7 M 0 and the minimum mass for 
neon ignition in the core is traditionally referred to as 
M mas « 10 M 0 (Becker & Iben 1979, 1980; Garcfa-Berro 
et al. 1997). Stars with zero age main-sequence (ZAMS) 
masses between ~ 7 M 0 and ~ 10 M 0 are designated as 
super-AGB stars (henceforth SAGB, Ritossa et al. 1996, 
1999; Gil-Pons et al. 2005; Siess 2006, 2007, 2010; Poe- 
larends et al. 2008; Doherty et al. 2010). Due to the 
inferred slope of the stellar initial mass function from 
observations (e.g., Jennings et al. 2012), single stars in 
this ZAMS mass range represent the population of stars 
that can produce the most massive white dwarfs, the 
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most numerous supernovae and possibly the least mas¬ 
sive neutron stars (e.g., Doherty et al. 2015). SAGB stars 
may also make significant contributions to the Galactic 
inventory of isotopes such as 7 Li, 14 N, 23 Na, 25_26 Mg, 
26_27 A1 and 60 Fe (Siess 2010; Ventura et al. 2013; Do¬ 
herty et al. 2014a,b). 

After helium is exhausted in the core, stars ascend¬ 
ing the SAGB develop partially electron degenerate car¬ 
bon-oxygen cores ranging from ^0.9 M 0 to ~ 2.0 M 0 
(pioneering studies of CO cores include Rakavy et al. 
1967; Beaudet & Salpeter 1969; Boozer et al. 1973). 
Depending primarily on the ZAMS mass but also sen¬ 
sitively on the composition mixing model (Poelarends 
et al. 2008; Siess 2009), the ignition of carbon may 
not occur at all (for stars < 7 M 0 ), occur at the cen¬ 
ter of the star (for stars > 10 M 0 ), or occur some¬ 
where off-center. In the off-center case, ignition is fol¬ 
lowed by the inward propagation of a subsonic burn¬ 
ing front (Nomoto & Iben 1985; Timmes et al. 1994; 
Garcfa-Berro et al. 1997; Saio & Nomoto 1998). Trail¬ 
ing behind the burning front is a convective region that 
may extend outward « 0.6 M 0 ; see Figure 1 for an il¬ 
lustration. When a steady-state, convectively bounded, 
subsonic, carbon burning front (henceforth a “flame”) 
propagates toward the center of the CO core, the flame 
leaves behind oxygen-neon-magnesium (ONeNa) ashes. 
The inward propagating carbon flame may or may not 
reach the center of the star, depending on the parame¬ 
ters adopted for composition mixing beyond the convec¬ 
tive boundary set by mixing-length theory (e.g., thermo- 
haline, overshoot, semi-convection, Siess 2009; Stancliffe 
et al. 2009; Denissenkov et al. 2013). If the flame makes 
it to the center, then the original CO core is converted 
into an ONeNa core. Such SAGB stars can explode as 
electron capture super novae if their ONeNa core masses 
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reach central densities in excess of the threshold density 
for the 20 Ne(e _ ,^) 20 F electron capture reaction (Miyaji 
et al. 1980; Nomoto 1984; Gutierrez et al. 1996; Jones 
et al. 2013), as may be the case for the Crab Nebula 
(Davidson et al. 1982; Nomoto et al. 1982; Wanajo et al. 
2009) or for potentially explaining observations of sub- 
luminous type II-P supernovae (Smartt 2009). If the 
flame does not make it to the center, then the star is 
left with inner parts of the original CO core surrounded 
by a layer of ONeNa. Such hybrid white dwarfs may 
provide unusual Type la supernovae progenitors (Siess 
2009; Denissenkov et al. 2013; Wang et al. 2014; Chen 
et al. 2014). 
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Fig. 1.— Illustration of the structure of an SAGB star during 
carbon burning. In the center is carbon/oxygen degenerate core, 
surrounded by a layer of helium which is then surrounded by a 
hydrogen envelope. Ignition of carbon is followed by a trailing 
convective region that drives a flash or flame towards the center. 
The scale on the left is the mass density, and the scale on the right 
is the enclosed mass. 

This paper explores the ignition and subsequent evolu¬ 
tion of carbon burning in SAGB stellar models as a func¬ 
tion of the ZAMS mass, initial rotation rate, the mag¬ 
nitude of various mixing parameters such as convective 
overshoot, semiconvection, thermohaline and angular 
momentum transport. We sample this multi-dimensional 
space with 2,755 MESA stellar evolution models (Pax¬ 
ton et al. 2011, 2013) that are evolved from the pre- 
main-sequence to the end of carbon burning. All models 
have Z = 0.02 and a solar composition from Grevesse & 
Sauval (1998). Along the way we provide quasi-analytic 
models for interpreting the results. In § 2 we discuss the 
input physics for our calculations, the composition mix¬ 
ing processes considered, and the MESA implementation 
of rotation and magnetic fields. In § 3 we define the base¬ 
line parameters used for our calculations and discuss the 
grids used for our exploration into this three dimensional 
parameter space. In § 4 we present the results of our non¬ 


rotating models, an analytical approximation of the loca¬ 
tion of first carbon ignition, and the evolution of carbon 
burning flames and flashes. In § 5 we present the results 
of the effect of rotation and overshoot on the ignition, 
evolution, and death of carbon burning. In § 6 we study 
the impact of the semiconvection, thermohaline, and an¬ 
gular momentum transport coefficients on the location 
of first carbon ignition on our results. In § 7 we present 
the results of spatial and temporal convergence studies 
on our results, and in § 8 we discuss our results and their 
implications. 

2. INSTRUMENT AND METHODS 

Our numerical instrument is MESA version 6794. We 
use the included sagb_NeNa_MgAl.net reaction network, 
which follows 22 isotopes from 1 H to 27 Al to track hy¬ 
drogen (pp chains, CNO-, NeNa-, and MgAl-cycles), he¬ 
lium and carbon burning. The 51 thermonuclear reac¬ 
tion rates coupling these isotopes are from JIN A reaclib 
version V2.0 2013-04-02 (Cyburt et al. 2010), energy- 
loss rates and their derivatives from thermal neutrinos 
are from the fitting formulae of Itoh et al. (1996), and 
electron screening factors for thermonuclear reactions in 
both the weak and strong regimes are from Dewitt et al. 
(1973); Graboske et al. (1973) and Alastuey & Jancovici 
(1978) with plasma parameters from Itoh et al. (1979). 
Poelarends et al. (2008) showed that increasing the mass 
loss rate could decrease the mass range for systems that 
will become electron capture supernovae (ECSNe). We 
thus use a Reimer mass loss prescription (Reimers 1975) 
with r]=0 .5 on the RGB and a Blocker mass loss pre¬ 
scription (Bloecker 1995) with 77 = 0 .05 on the AGB. The 
MESA inlists are publicly available 3 . 

Analysis of a carbon burning event requires knowledge 
of when and where carbon ignites. We identify the cell lo¬ 
cation of carbon burning, hence the ignition mass coordi¬ 
nate (M/ jS ), by three criteria. First, we require e nuc e v 
in a CO core. Second, we require that 4 He is depleted 
in the ignition region as some of the lowest mass stars 
investigated would have a small amount of 12 C burning 
near the CO core, 4 He shell boundary. Finally, we re¬ 
quire X( 20 Ne) > X( 23 Na) > X( 24 Mg), which indicates 
we have vigorous 12 C -f 12 C burning. 

The end of a carbon burning event is defined when no 
cell within 10% of the mass location of the flame, dur¬ 
ing the next time step, has e nuc ^ e y . We then define 
the final flame location (M/ ?e ), where carbon burning 
is quenched, as the minimum value of the mass loca¬ 
tion taken at the end of all the carbon burning events. 
This is independent of any subsequent carbon flashes, as 
we record only the closest approach the carbon burning 
makes to the core. 


2.1. Mixing 

Treatment of convective processes within stellar inte¬ 
riors is essential for a physically accurate stellar model. 
We briefly discuss the composition mixing processes used 
in our calculations, how MESA models the mixing pro¬ 
cesses, and previous studies that guide our choices for 
our baseline mixing parameters. Values for our baseline 
parameters are given in § 3. 

3 http://mesastar.org/results 
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We use the Schwarzschild criterion for convection along 
with the Cox implementation of Mixing Length Theory 
(MLT) (Cox & Giuli 1968). The Schwarzschild criterion 
describes that a region is stable to convection if the gra¬ 
dient of a piece of adiabatic matter is less than that of 
the temperature gradient of the stellar atmosphere: 


din T 
d In P 


> 

ad 


d In T 


d In P 


rad 


(i) 


MLT has a free parameter, ctmlt, as described by Bohm- 
Vitense (1958). Values within 1.6 < c^mlt 2.2 have 
been inferred by Noels et al. (1991) who compared obser¬ 
vations of the a Centauri binary star system with stellar 
evolution models. Trampedach et al. (2014) also suggests 
1.6 < c^mlt < 2.05 from calibrating ID stellar models to 
3D radiation-coupled hydrodynamics simulations of con¬ 
vection in stellar surface layers. 

Turbulent velocity fields have been suggested to de¬ 
cay exponentially beyond the Schwarzschild convective 
boundary defined by equation 1 (Herwig et al. 1997; Ven¬ 
tura et al. 1998; Mazzitelli et al. 1999; Blocker et al. 2000; 
Herwig 2000) leading to a diffusive treatment of mixing 
as 

D 0 v = D con ^oexpf-— T—) . (2) 

V Jov^PflJ 


Here D conVj o is the convective diffusion coefficient at 
the convective boundary, z is the radial distance from 
the convective boundary, Ap ? o is the local pressure scale 
height and / ov is an adjustable parameter. MESA offers 
the flexibility of allowing / ov to be different for different 
convective regions (H burning, He burning, metal burn¬ 
ing and non burning). However, we set / ov to the same 
value in all convective regions (see §3 for the values cho¬ 
sen). 

Semiconvection occurs when regions of the stellar inte¬ 
rior are stable to the Ledoux criterion and unstable to the 
Schwarzschild criterion (Kippenhahn et al. 2012). This 
occurs when V a ^ < Vp < Vl, where Vl = V a ^ + B 
and B is the Brunt composition gradient. MESA treats 
semiconvection as a diffusive process (Langer et al. 1983, 
1985; Heger et al. 2000; Zaussinger & Spruit 2013) with 
a diffusion coefficient 


A* — &SC 
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where K is the radiative conductivity, C v is the specific 
heat at constant pressure and a sc is an adjustable di¬ 
mensionless parameter describing the speed with which 
convective mixing may occur at the boundary defined by 
equation 1. Ongoing efforts to calibrate such semiconvec¬ 
tion models include multidimensional numerical simula¬ 
tions of double-diffusive convection (Zaussinger & Spruit 
2013; Spruit 2013) and comparing massive star models 
with observations (Yoon et al. 2006). 

Thermohaline mixing occurs when Vt — V a d < B < 0. 
These are regions stable against convection, according to 
the Ledoux criterion, but have an inversion of the mean 
molecular weight (Kippenhahn et al. 1980). This type 
mixing forms elongated fluid parcels, sometimes called 
“salt-fingers”. MESA treats thermohaline mixing as a 
diffusion process (Ulrich 1972; Kippenhahn et al. 1980; 


Brown et al. 2013; Zemskova et al. 2014) 

_ 3 KB 
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where a t h is a dimensionless parameter, related to the 
aspect ratio of the salt fingers. Estimates for this pa¬ 
rameter range from 1 < a t h < 667 (Kippenhahn et al. 
1980; Charbonnel & Zahn 2007; Cantiello & Langer 2010; 
Stancliffe 2010; Wachlin et al. 2011), with some multidi¬ 
mensional simulations suggesting this parameter is signif¬ 
icantly overestimated (Denissenkov 2010; Traxler et al. 
2011; Denissenkov & Merryfield 2011; Lattanzio et al. 
2015). 


2.2. Rotation and Magnetic Fields 

MESA implements rotation by making the assumption 
that the angular velocity, ca, is constant over isobars; see 
Paxton et al. (2013) for the implementation of rotation 
into MESA. Such an assumption is often referred to as 
shellular approximation (Zahn 1992; Meynet & Maeder 
1997), and allows the stellar structure equations to be 
solved in one dimension for a rotating star. For this 
study, rotation is initialized by imposing a solid body 
rotation law at ZAMS, where the total luminosity equals 
the nuclear burning luminosity. 

The transport of angular momentum and material 
due to rotationally induced instabilities is followed us¬ 
ing a diffusion approximation (e.g., Endal & Sofia 1978; 
Pinsonneault et al. 1989; Heger et al. 2000; Maeder 
& Meynet 2003, 2004; Heger et al. 2005; Suijs et al. 
2008) for the dynamical shear instability (DSI), secu¬ 
lar shear instability (SSI), Eddington-Sweet circulation 
(ES), Goldreich-Schubert-Fricke instability (GSF), and 
Spruit-Tayler dynamo (ST). See Heger et al. (2000) for a 
description of the physics of the different instabilities and 
the calculation of the respective diffusion coefficients. 

Berger et al. (2005) investigated the Ca line profiles of 
a sample of DA white dwarfs, concluding their rotational 
velocities are generally less than 10 km s -1 . These values 
are significantly less than the values determined by the 
rotating, nonmagnetic models of Langer et al. (1999). 
Internal magnetic torques as proposed by Spruit (1998) 
have been suggested as an effective mechanism to spin 
down the cores of these white dwarf progenitors during 
the giant phase. Suijs et al. (2008) showed that magnetic 
torques as calculated in Spruit (2002) produce rotational 
velocities in better agreement with the observed values 
of Berger et al. (2005). We therefore include internal 
magnetic fields and the Spruit Tayler dynamo angular 
momentum mechanism for our rotating MESA models. 

Magnetic fields are implemented in MESA using the 
formalism of Heger et al. (2005), where a magnetic torque 
due to a dynamo (Spruit 2002) allows angular momen¬ 
tum to be transported inside the star. The radial com¬ 
ponent, B r , and the azimuthal component, of the 
magnetic field are modeled as 

Bj, ~ r(4wp) 1/2 u>A (5) 

B r ~ B^rky 1 , (6) 

where r is the radial coordinate, p the density, oj_\ the 
Alfven frequency, and k the wavenumber. These mag- 
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Fig. 2. — Three slices explored in the (mass, rotation, overshoot) parameter space. We calculate 1326 models in mass-rotation rate plane 
(blue slice), 546 models in the rotation rate-overshoot plane (red slice), and 546 models in the mass-overshoot plane (green slice). 


netic fields then provide a torque 


which acts to slow down the stars rotation rate by de¬ 
creasing the amount of differential rotation inside the 
star (Heger et al. 2005). 

The initial rotation is normalized against the critical 
rotation rate for the star Q C rit = y/(l — T/ L e dd)cM / R 3 , 
where c the speed of light, M the mass of the star, R the 
stellar radius, L the luminosity and L e dd the Eddington 
luminosity. The initial magnetic field, B r = B$ = 0 for 
all our rotating models. Effects of rotationally induced 
mass loss are not included. 

3. GRIDS 

We define a set of baseline parameters and construct 
a number of grids surrounding that baseline set to in¬ 
vestigate the sensitivity of carbon burning in the SAGB 
models with respect to variation in the baseline parame¬ 
ters. Choices in the numerical values of the baseline pa¬ 
rameters are based on the current understanding of the 
canonical values for SAGB and other stars when using 
MESA. Choices in the range of values explored were de¬ 
signed such that we could explore the various competing 
factors involved in stellar evolution models. 

Table 1 lists the baseline mixing, spatial resolution, 
and temporal resolution parameters. The parameter f c is 
the ratio of turbulent viscosity to the diffusion coefficient 
(Heger et al. 2000), f^ is the sensitivity to /r-gradients 
(Heger et al. 2000), while £mesh controls the spatial reso¬ 
lution by determining the relative magnitude of changes 
between the adjacent cells (Paxton et al. 2011), and w t 
controls the temporal resolution by modulating the mag¬ 
nitude of the allowed changes between time steps (Paxton 
et al. 2011). Baseline values for ^ mes h and w t were based 


on computational requirements. However see § 7 for a 
discussion into their relative effects on our results. The 
other baseline parameters listed in Table 1 are discussed 
in §2.1. 

3.1. Mass-Rotation-Overshoot Grid 

We explore the (ZAMS mass, initial rotation, over¬ 
shoot) parameter space with three slices through this 3D 
data cube. Table 2 lists the start, stop, and step values 
for two of the three quantities while holding the third 
quantity constant. The number of SAGB models is 1326 
in the mass-rotation plane, 546 in the rotation-overshoot 
plane, and 546 in the mass-overshoot plane for a total 
of 2418 SAGB models. Figure 2 illustrates these three 
orthogonal slices though this 3D parameter space. 

For each slice within the 3D data cube we chose one 
quantity to be held fixed. In the mass-rotation plane 
this is the overshoot value, and our choice of f ov =0.016 
is based on the canonical value for this overshoot model 
(Herwig 2000). In the rotation-overshoot plane we held 
the initial mass fixed at M = 8M 0 , which was se¬ 
lected because, based on the non-rotating models, we 
expected we could induce a range of behaviors from non¬ 
ignition, to off-center ignition to central ignition. In 
the mass-overshoot plane, we held the initial rotation 
fixed at (fl/Q cr ^)i=0.25, purely as a middle ground value 
between non-rotating models and our upper bound of 
(fl/fl cr ^)i=0.5. Our choice for the range of values cov¬ 
ered by each grid was based on a requirement to have a 
comprehensive sample over the canonical baseline values 
for SAGB stars. 

3.2. Mixing Coefficients Grid 

We investigate the (semiconvection, thermohaline, 
overshoot, angular momentum diffusion) parameter 
space on the location of the first carbon ignition in 
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TABLE 1 

Baseline Parameters 


Parameter 

Value 

Mixing Length Theory (ckmlt) 

2.0000 

Semiconvection (a sc ) 

0.0100 

Thermohaline (a t h) 

1.0000 

Overshoot (/ov) 

0.0160 

Angular Momentum ( r] a m ) 

1.0000 

Turbulent Viscosity (/ c ) 

0.0333 

/^-gradient Sensitivity (/^) 

0.0500 

Mesh Delta Coefficient (5 m esh) 

0.5000 

Variance Control Target (wt) 

0.0001 


TABLE 2 

Mass-Rotation-Overshoot Grid 


Variable 

Start 

Stop 

Step 

Constant 

Mzams 

6.0 

11.0 

0.1 

/ou =0.016 

(D / Tl C rit)i 

0.0 

0.5 

0.02 


(Q / Q cr it)i 

0.0 

0.5 

0.02 

^zams=8.0 

fov 

0.0 

0.02 

0.001 


Mzams 

6.0 

11.0 

0.2 

( Q / n crit ) i = 0.25 

fov 

0.0 

0.02 

0.001 



TABLE 3 

Mixing Coefficients Grid 


Variable 


Values 


CX-sc 

0.0 

10“ 3 

10“ 2 

10- 1 

<*th 

0.0 

0.1 

1.0 

10.0 

fov 

0.000 

0.001 

0.016 

0.020 

Vam 

0.0 

0.5 

1.0 

1.5 


TABLE 4 

Spatial and Temporal Convergence Grid 

Variable 


Values 


^ZAMS 

7 

8 

9 

{Q / Q cr it)i 

0.0 

0.25 

0.5 

^mesh 

0.1 

0.5 

1.0 

w t 

10- 5 

l 

0 

T—1 

10“ 3 


a 8 M 0 ZAMS, (Q/Tl cr i t )i=0.2b model with selected 
points in this 4D data cube. Table 3 lists these quan¬ 
tities and their selected values. We choose to limit the 
range of a t h \n this grid to span only the lower values dis¬ 
cussed in the literature. rj arn is a scale factor that alters 
the strength of angular momentum diffusion in each cell; 
see Paxton et al. (2013) for details. The total number 
of SAGB models in this grid is 4 4 =256. This grid per¬ 
mits assessment of the relative strengths of each mixing 
process. 

3.3. Spatial and Temporal Convergence Grid 

Finally, we examine the spatial and temporal conver¬ 
gence properties of a subset of our SAGB models. The 
MESA parameter ^ mes h broadly controls the spatial reso¬ 
lution and ^broadly relates to the temporal resolution. 
For £ m esh=0.5 there are « 5,000—10,000 spatial points 
center to surface. Spatial resolutions necessary to cap¬ 


ture carbon burning flames and flashes are discussed in 
§7. For w t = 10 -4 a temporal resolution of ~ 10 yrs 
is achieved. Table 4 lists these quantities and their val¬ 
ues. The total number of SAGB models in this resolution 
sensitivity grid is 3 4 =81. 

4. RESULTS FROM NON-ROTATING MODELS WITH 
BASELINE MIXING PARAMETERS 

We begin by considering a series of non-rotating stel¬ 
lar models using our baseline mixing parameters as de¬ 
scribed in § 2. Figure 3 shows the Kippenhahn plots 
of these non-rotating stars as representative samples for 
all the stellar models. In the 7 M 0 case, carbon ignites 
off-center at ~ 0.6 Mq. However, it rapidly quenches 
and does not propagate towards the core. The 7.5 M 0 
case undergoes a series of carbon flashes, with each flash 
igniting closer to the core but it is unable to form a 
steady state flame. In the 8 M 0 case, an off-center igni¬ 
tion occurs at ~ 0.15 M 0 . A flame propagates inwards 
and through a series of distinct flashes as the flame ap¬ 
proaches the core, and almost reaches the center. For 
the 9 M 0 case, carbon ignites at the center. In both the 
8 and 9 M 0 models, secondary flashes at ~ 0.5 M 0 are 
due to off-center carbon burning. For all cases the core is 
undergoing significant cooling, primarily through photo¬ 
neutrino and plasma neutrino losses prior to the first ig¬ 
nition of carbon (Nomoto 1984, 1987; Ritossa et al. 1996, 
1999). The 7.5, 8.0 and 9.0 M 0 cases the stars undergo a 
series of subsequent carbon flashes that travel away from 
the core. For the 9.0M 0 case with no overshoot (bottom 
right in Figure 3) the flame ignites off-center, contrary to 
the center ignition of the 9 M 0 model with / ov . The star 
undergoes a flash and then a flame, which propagates all 
the way to the center. The model with closest morphol¬ 
ogy is the 8 Mq case which has a flash and flame, but 
carbon burning does not reach the center. The differ¬ 
ence between models with overshoot and without will be 
discussed further in §6. 

The He and CO core masses for the non-rotating mod¬ 
els are shown in Figure 4 as a function of the ZAMS mass. 
Stellar models with Mzams < 7 Af 0 do not ignite carbon 
(Becker & Iben 1979, 1980) and are not shown. The car¬ 
bon core mass increases linearly with ZAMS mass over 
the range shown. For Mzams ^ 7.5 M 0 , the He forms 
a radiative, geometrically thin, burning shell accreting 
material onto the CO core at a rate of ss 1O -6 M 0 yr -1 . 
Between ss 7.5 M 0 and ~ 8.0 M 0 , Figure 4 shows the 
He shell transitions from an geometrically thin shell to an 
geometrically thick shell. For Mzams <8.0 M 0 , the geo¬ 
metrically thick He shell mass grows linearly with ZAMS 
mass over the range shown, accreting material onto the 
CO core at ss 3xlO -6 M 0 yr -1 at 8.0 M 0 , and increasing 
to ss 2.OxlO -5 M 0 yr -1 for the 11 Mq stellar models. 

As the star evolves, the He shell grows in mass reach¬ 
ing ^1.7 Mq for the 7 Mq stars and up to 3.2 for the 
11 Mq stars. Once the He shell reaches its peak size, 
the 4 He depletes quickly leaving a CO core. Shortly af¬ 
terwards the 4 He shell begins rapidly accreting material 
onto the CO core, eventually depleting itself as well. At 
the lowest ZAMS masses (< 7M 0 ), ignition occurs af¬ 
ter the 4 He shell has accreted onto the CO core and the 
star has finished its second dredge up (2DU) (Becker & 
Iben 1979). In the transition region where we have a se- 
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Fig. 3.— Kippenhahn diagrams of flashes and steady state flames of 7.0 M© (top left), 7.5 M© (top right), 8.0 M© (middle left), 8.2 M© 
(middle right) and 9.0 M© (bottom left) ZAMS model stars, for non-rotating f ov = 0.016. Bottom right a 9.0 M© non-rotating, / ov =0.0 
model. Dark purple regions indicate regions of cooling, primarily from thermal neutrino losses, with the darker shades of purple representing 
a logarithmic increase in the cooling rate. Red regions indicate significant nuclear burning, light blue regions indicate convection. For 
clarity the regions undergoing other types of mixing are not shown. 
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ries of flashes (~ 7 — 7.8M 0 ), the stars are undergoing 
their 2DU while igniting carbon. In the higher mass sys¬ 
tems, « 7.9—8.2 M 0 , where we have steady state flames 
or central carbon ignition, the star ignites carbon before 
the 2DU and before any significant accretion on to the 
CO core can occur. Above 8.3 M 0 we have dredge out 
events (Ritossa et al. 1999), where the 4 He shell grows 
an an outwardly moving convection zone which merges 
with the inwardly moving convective envelope. 

At the base of the flame, we can ask whether a packet of 
convective material can penetrate into the region ahead 
of the flame transferring energy which will decrease the 
flames lifetime. A simplified derivation (L. Bildsten, pri¬ 
vate communication) assumes the length scale l of the 
flame front (the distance over which the temperature de¬ 
creases from a peak inside the flame to the background 
value) is much less than the local pressure scale height, 
l <C H. This implies the pressure is constant across the 
subsonic flame front. A convective packet will move from 
a region of high temperature to a lower temperature re¬ 
gion, at constant pressure. Assuming adiabatic motion, 
a fluid packet keeps its original temperature. Thus, the 
buoyancy felt by this convective packet is: 


a = 


dlnT 

g——Sr . 
dr 


(8) 


Where a is the buoyancy acceleration, g is the local accel¬ 
eration due to gravity, and Sr is the distance the packet 
moves ahead of the flame. Simplifying, 


a 


—i sr • 


( 9 ) 


and solving this harmonic oscillator equation we find 


Sr ~ 


Vc 

N ’ 


( 10 ) 


where v c is the convective velocity and N is the 
Brunt-Vaisala frequency. Using the local scale height 
H = c 2 s /g, where c s is the local sound speed: 


Sr _ v c /iJ\ 1/2 

T “ ^ vTy 


(ii) 


In a typical carbon flame we have v c /c s « lO -4 and 
H/l « 10, thus Sr <C l and a convective fluid packet 
cannot penetrate the flame front. 


4.1. First Ignition of Carbon 

Figure 5 shows the mass location of the first ignition 
as a function of the stellar model ZAMS mass. For our 
choice of baseline mixing parameters, the lowest mass 
star to ignite carbon is a 7 M 0 model and ignition oc¬ 
curs off-center at a mass coordinate of « 0.65 M 0 . As the 
ZAMS mass increases, the location of the off-center first 
ignition location moves steadily inwards in mass (Siess 
2007). For ZAMS masses larger than 8.4 M 0 , carbon ig¬ 
nites at the center. Figure 5 also shows the local density 
at the location of first ignition. All stellar models that ig¬ 
nite carbon off-center, 7 M 0 < Mzams < 8.4 M 0 , do so at 
a nearly constant density of logp ~ 6.2, or p « 1.5xl0 6 
g cm -3 . For stellar models that undergo central car¬ 
bon ignition, 8.4 M 0 < Mzams <11 Mq, the density 
at ignition monotonically decreases. Models which do 
not ignite carbon will eventually form a CO WD. Those 



Fig. 4. — He and CO core mass as a function of ZAMS mass at 
the first ignition of carbon , for our baseline mixing parameters. In 
blue is the 4 He core mass, while in red is the CO core mass. Stars 
with Mzams < 7.0 M 0 are not shown as they do not have a carbon 
ignition point. SAGB models with 7 M 0 < Mzams < 7.5 M 0 have 
thin helium envelopes while models with Mzams >7.5 M 0 have 
thick helium envelopes. 



Mzams [M©] 


Fig. 5. — Mass location of first carbon ignition (red) and local 
density at first carbon ignition (blue) as a function of the ZAMS 
mass for non-rotating, baseline mixing parameters models. 

models that undergo off-center carbon ignition but where 
the burning does not reach the center will form hybrid 
CO+ONe WD. Model stars which ignite carbon at the 
center will eventually form an ONeNa WD which may 
explode as an ECSNe. 

To a first approximation, at ignition the nuclear burn¬ 
ing timescale and thermal diffusion timescales are equal 
(Timmes & Woosley 1992) 

a E 

Tdiff ~ —7T T burn ~ - > (12) 

P C p ^nuc 

where a is the thermal conductivity, C p is the specific 
heat capacity at constant pressure, E is the local thermal 
energy, and e nuc is the screened nuclear energy generation 
rate during carbon burning including energy losses due 
to neutrino cooling. For a given temperature, density, 


log Pf,s [gem 3 ] 








and composition both C v and E may be calculated from 
an equation of state. For carbon burning, e nuc takes the 
form (e.g., Woosley et al. 2004): 


6.7 X ICt 3 X 2 ( 12 C) P6 /screen Ai 2> 12 


(13) 


where X( 12 C) is the carbon mass fraction, p 6 is the den¬ 
sity divided by 10 6 g cm -3 , \ 12 p 2 is the unscreened nu¬ 
clear reaction rate for carbon burning and /screen is the 
electron screening factor. Using a MESA 501 isotope re¬ 
action network that includes neutrino losses with an ini¬ 
tial initial composition of X( 12 C) = 0.3 and X( 16 0) = 0.7 
to calculate e nuc over the relevant range in the p -T plane, 
we find positive values of e nuc may be approximated by 
the power law 


1.6 x 10 7 


7 x 10 8 


23 


2 x 10 6 


1.2 


(14) 


Results for carbon ignition for any large reaction net¬ 
work, including the 501 isotope network used here, are 
expected to be similar to that of the smaller 22 isotope 
net used in the SAGB models because both networks 
have the same set of key isotopes and reaction rates cru¬ 
cial for carbon burning, (e.g., Timmes et al. 2000). 
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Fig. 6.— Nuclear burning (top) and thermal diffusion (bottom) 
timescale in the p-T plane for X( 12 C) = 0.3 and X( 16 Q) = 0.7. 


Fitting the other quantities in equation 12 in a similar 
manner, we find the following expressions for the nuclear 
burning timescale and thermal diffusion timescales 


Tburn = 5.1 X IQ 9 


T dif f = 4.0 X 10 9 


7 x 10 8 


7 x 10 8 


-32 


- 2.4 


P 


2 x 10 6 


P 


2 x 10 6 


- 0.8 


-1 


(15) 


These two timescales are shown in Figure 6. 
the two timescales gives 


7 x 10 8 


29.6 


2 x 10 6 


- 0.2 


= 1.3 


• (16) 
Equating 


(17) 


At the threshold of vigorous carbon burning, 
T « 7xl0 8 K, this expression gives a unique and 
constant ignition density 


Pigr 


2.1 x Kr g cm 


-3 


(18) 


which is consistent with the constant ignition density of 
p « 1.5 xlO 6 g cm -3 found in the MESA SAGB models 
that ignite carbon off-center, 7 M 0 < Mzams < 8.4 M 0 , 
of Figure 5. We find this result also holds for our rotating 
SAGB models. 



Fig. 7.— Polytrope fits to the MESA carbon core density struc¬ 
ture for the non-rotating (a) 7.0 MQ,(b) 8.0 Mq, and (c) 9.0 Mq 
ZAMS model stars from top to bottom, respectively. In each case, 
the MESA model is bounded by two polytropic fits; (a) red n=2.9, 
blue n=3.1; (b) red n=3.6, blue n=3.8; (c) red n=4.2, blue n=4.4. 

Figure 7 (a) shows the carbon core of a non-rotating 
7 M 0 ZAMS star can be well approximated by a n « 3.0 
polytropic model, (b) a 8 M 0 ZAMS star by a n « 3.7 
polytrope, and (c) a 9 M 0 ZAMS star by a n ~ 4.3 poly¬ 
trope. The density structures for these polytropic models 
were calculated using an open-source tool 4 . In addition, 
the polytropic relations offer an explanation for why the 
location of the first ignition moves steadily inwards in 

4 http://cococubed.asu.edu/code_pages/polytrope.shtml 
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mass for the off-center ignition cases (see Figure 5). The 
density structure of the n = 3.0, 3.7, and 4.3 polytrope 
models are shown in Figure 8. The mass locations for 
a fixed ignition density (dashed horizontal line), where 
ignition occurs (dash vertical lines) moves monotonically 
inwards as the polytropic index increases, as the carbon 
core mass increases, as the ZAMS mass increases. 



Fig. 8.— Mass locations of first carbon ignition from the ignition 
curve and polytropic models in the density—mass plane. 


We now turn to the decrease in the central density for 
those SAGB models in Figure 4 that centrally ignite car¬ 
bon. Homology relations between the central density p c 
and the mass M (e.g., Hansen et al. 2004; Kippenhahn 
et al. 2012) for a chemically homogeneous star charac¬ 
terized by a mean molecular weight /i, constant opacity, 
ideal gas equation of state and power-law energy gener¬ 
ation rate e nuc oc T v indicate 

3(4—i/) _ _ 2(3— v>) 

p c oc p ^+ 3 M "+ 3 . (19) 

For carbon burning near ignition, equation 14 shows 
v « 23, and equation 19 then gives p c oc M -1,5 . The 
negative exponent shows that the density at first igni¬ 
tion, for those 8.4 M 0 < Mzams <11 M 0 models that 
undergo central carbon ignition, monotonically decreases 
as the mass of the carbon core increases with a slope that 
is consistent with the rate of decline shown in Figure 5. 

4.2. Carbon Burning flashes and Transition to a Steady 

State Flame 

Local characteristics before, during, and after the first 
off-center carbon flash in a 7.5 M 0 model is shown Fig¬ 
ure 9. Before the first ignition of carbon, the CO core 
is in its most compact, most electron degenerate config¬ 
uration. The first off-center carbon burning flash is thus 
the most powerful; any subsequent flashes or steady-state 
flames take place under more extended, less degenerate 
configurations. In addition, the energy released during 
the first ignition decreases as the ZAMS mass increases 
because the CO core is not as compact and not as degen¬ 
erate. Furthermore, the first carbon flash impacts the 
base of the convective envelope more strongly in lower 
mass SAGB models than in higher mass SAGB models 
(Garcfa-Berro et al. 1997; Siess 2006), partly because of 


their more compact configuration but also because the 
first flash occurs farther from the center in lower mass 
models (see Figure 5). 



Fig. 9. — Local and global characteristics before, during, and 
after the first carbon burning flash in an 7.5 Mq model. L 0 is the 
total luminosity due to carbon burning, L# e the luminosity due to 
nuclear helium burning, L u the luminosity due to thermal neutrino 
losses, T c the central temperature, p c central density, R the surface 
radius, and L s the surface luminosity. The black dotted line marks 
first ignition of carbon. 

At the start of carbon burning at the ignition density 
of Pign ~ 1.5xl0 6 g cm -3 , marked by the vertical dotted 
line in Figure 9, the energy generation rate rapidly rises. 
Figure 9 shows that « 10% of the energy produced by 
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nuclear reactions is lost to neutrinos, with the balance of 
the thermal energy transported by convection (Ritossa 
et al. 1996; Saio & Nomoto 1998). This is not the « 50% 
expected for steady state burning, as the 7.5M 0 model 
is undergoing a time-dependent flash and not a steady 
state flame. Were it not for carbon burning, the sur¬ 
face luminosity and radius of these SAGB models would 
continue to increase in step with an increasing helium- 
burning luminosity, similar to lower mass, Mzams 7 
M© AGB stars that do not experience carbon burning 
(Garcfa-Berro et al. 1997; Karakas & Lattanzio 2014). 
These dramatic events mainly impact the core and have 
only a modest effect on the structure of the outer enve¬ 
lope and surface luminosity. 

Figure 9 shows the rapid injection of energy causes 
the entire core to expand and cool (Ritossa et al. 1999). 
A convective region also develops above the ignition lo¬ 
cation in response to the injection of energy from car¬ 
bon burning. The overall expansion of the core decreases 
the electron degeneracy parameter (rj = p e /kT — oo 
means non-degenerate where fi e is the electron chemi¬ 
cal potential, r] oo means perfect degeneracy). The 
expansion extinguishes the flash partially because the 
nuclear energy generation rate drops below the critical 
luminosity (e.g., Siess 2006) and partially because the 
ignition density pi gn is pushed deeper into the stellar 
model. The first flash thus quenches and does not be¬ 
come a steady-state flame that propagates inwards to¬ 
ward the center. 

After the first flash is quenched and the deposition of 
energy from nuclear reactions ceases, the core again con¬ 
tracts but to a less compact, less degenerate configura¬ 
tion. This contraction leads to the second ignition of 
carbon. The second flash (and any subsequent flashes) 
occurs at a deeper mass coordinate because the accre¬ 
tion of C from the He burning shell slightly increases 
the core mass which moves the location of the critical 
density pi gn inwards where there is also fresh, previously 
unburned fuel. Similar evolution pathways can occur for 
other fuels and other masses, for example, neon burning 
in more massive models (e.g., Jones et al. 2013). 

The evolution of the second and any subsequent flashes 
is notably different than the first flash as shown in Fig¬ 
ure 3. When these later flashes develop, their convective 
region grows into regions previously occupied by previ¬ 
ous flashes where carbon had been depleted. Thus the 
nuclear energy production rate, which from equation 13 
is proportional to X 2 ( 12 C), is reduced. Furthermore, the 
subsequent core expansion and induced thermodynamic 
changes are significantly smaller. The first flash lasts the 
longest and is the most vigorous, while subsequent flashes 
have shorter lifespans, release less energy, and the time 
interval between flashes becomes shorter Models which 
have a geometrically thin helium envelope show more 
flashes than models with a geometrically thick helium 
envelope (see Figure 4). 

Each flash releases less energy, expands the core by a 
smaller amount, and moves the mass location of the crit¬ 
ical density inward at a slower rate. This allows a flash 
to transition to a steady-state flame. Thus, after one or 
more flashes a steady-state may be achieved. Combus¬ 
tion in steady-state flames is incomplete; only a small 
portion of the carbon burns. A condition of balanced 
power is set up where the rate of energy emitted as neu¬ 


trinos from the base of the convective region equals the 
power available from the unburned fuel that crosses the 
flame front. The inward propagation of the flame by ther¬ 
mal conduction is limited by the temperature at the base 
of the convective shell, which cannot greatly exceed the 
adiabatic value. These two local conditions give a unique 
speed for the flame, with typical values of ~ 0.1 cm s -1 . 
We verified the flame speed at several locations in the 
SAGB models and it is commensurate with previous lo¬ 
cal studies, with speeds of order 10 -3 — 10 _2 cm s -1 (e.g., 
Timmes et al. 1994; Ritossa et al. 1996; Garcfa-Berro 
et al. 1997; Siess 2006; Denissenkov et al. 2013). The 
flame lives ~ 20,000 years on its journey towards the 
center. 

The convective nature of the material behind the flame 
has two key consequences for its journey to the core. 
First, the temperature behind the flame front (i.e., to¬ 
wards the surface) is bounded. The ONeNa ashes of the 
burning are not allowed to assume an arbitrary value 
of temperature; rather convection fixes the temperature 
behind the flame front. The second feature of the convec¬ 
tive material is the lack of abundance gradients behind 
the flame. That is, convection uniformly mixes the ashes 
of the partial burning. 

For our standard mixing parameter settings, models 
with ZAMS masses in the 7—7.6 M 0 range do not achieve 
a steady-state flame. Instead, they undergo a series of 
flashes, where each flash occurs closer to the core (Siess 
2009; Denissenkov et al. 2013, and see Figure 3). The 
number of flashes increases as the ZAMS mass increases, 
until the ZAMS mass exceeds 7.6 M 0 when the first 
flash transitions into steady-state flame. For stars be¬ 
tween 7.7—8.4 M 0 the off-center steady-state flame be¬ 
gins closer to the center. The dependence of this mass 
range on the composition mixing parameters is discussed 
in §5. 

4.3. Does the Burning Reach the Center? 

Whether off-center carbon burning, either as a steady- 
state flame or as a series of time-dependent flashes, 
reaches the center in these models depends on the ZAMS 
mass and the adopted mixing parameters. When con¬ 
vective mixing operates within the Schwarzschild bound¬ 
aries, a flame will propagate all the way toward the center 
(Nomoto & Iben 1985; Garcfa-Berro et al. 1997; Saio & 
Nomoto 1998; Siess 2009; Doherty et al. 2010). Convec¬ 
tive overshoot is not strictly required for carbon flames 
to be quenched away from the center. Additional mixing 
processes such as thermohaline mixing have been shown 
to effectively quench the carbon flame away from the cen¬ 
ter while convection is operating within only the strict 
Schwarzschild boundaries. (Siess 2009; Denissenkov et al. 
2013) Mixing processes that, by design, extend beyond 
the MLT convective boundary take unburned carbon fuel 
ahead of the flame front and mix this fresh fuel with the 
ashes of the convective region behind the flame front. 
This starves the flame of fuel, with the nuclear energy 
production rate proportional to the square of the carbon 
abundance (see equation 13). For instance, instead of a 
fresh carbon mass fraction of 0.3 the carbon mass frac¬ 
tion near the ignition point may be depleted to 0.1 and 
polluted with enhanced abundances of 16 0, 20 Ne, 23 Na 
and 24 Mg. Examples of such convective boundary mix¬ 
ing processes include thermohaline and overshoot. For 
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large enough mixing parameters, the flame either disin¬ 
tegrates and sputters in a series of fuel starved flashes 
moving towards the center, or is extinguished before 
reaching the center (Siess 2009; Denissenkov et al. 2013). 
Where the flame is extinguished, if it is extinguished, can 
have repercussions for the composition of the subsequent 
white dwarf that is formed, and from there, possible con¬ 
sequences for supernova Type la models. 



Fig. 10. — Carbon burning extinction location in mass coor¬ 
dinate (red curve) and mass traversed (blue curve) for the non¬ 
rotating models with baseline mixing parameters. Stars with 
M > 9 Mq are not shown as they all ignite at the center. 

Figure 10 shows the location where carbon burning is 
extinguished (red curve, left y-axis) for the non-rotating 
models with baseline mixing parameters. The distance 
(in mass coordinates) that the flame traveled from birth 
to death is shown by the blue curve and the right y-axis. 
For 7.0 < Mzams < 8.2, carbon burning does not reach 
the center. As the ZAMS mass increases, the flame or 
flashes get closer to the center, eventually reaching the 
center at 8.2 M©. In terms of the mass traversed, the 
flame (or flashes) increases its travel distance from 7.0 
to 7.5 M© and then decreases for higher mass models. 
This transition occurs as stars with masses between 7.0 
and 7.5 M© undergo a one or more of flashes, where each 
flash does not travel, but each subsequent flashes ignites 
closer to the core. For higher mass models a steady-state 
flame reaches the center. The dependence of these results 
on the chosen mixing parameters is discussed § 5. 

For the models between 7.8 M© and 8.2 M©, which 
do not reach the center, they transition from a steady 
state flame into a series of flashes, as seen in Figure 3 
for the 8 M© case. These flashes ignite in regions where 
the X( 12 C) abundance has dropped due to the mixing. 
A subsequent flash thus requires a higher critical density 
(see equation 18) to ignite in the presence of a lower 
abundance of carbon. These flashes are able to drive a 
local expansion of the core and the critical density to 
move inwards into the core, analogous to the flashes seen 
in the < 7.8 M© models. The flashes eventually stop 
forming once the core can no longer reach the critical 
density leaving unburned carbon in the core. 



Fig. 11. — Angular momentum diffusion coefficients (top panel); 
specific angular momentum j, angular frequency uo (middle panel), 
and radial magnetic field B r , and azimuthal magnetic field B <© 
components (bottom panel), of an 8.0M© ZAMS model with an 
initial rotation at ZAMS of (£l/£} cr it)i= 0.2 , baseline mixing pa¬ 
rameters, at the onset of carbon ignition. The blue shaded region 
indicates a convective region, and the dashed (black) line shows the 
boundary of the CO core. The angular momentum diffusion coeffi¬ 
cients shown are the Eddington-Sweet circulation (ES), Goldreich- 
Schubert-Fricke instability (GSF), and Spruit-Tayler dynamo (ST). 


5. RESULTS FROM THE GRID OF ROTATING AND 
OVERSHOOTING MODELS 

In this section we investigate the impact of rotation 
and convective overshoot on carbon burning in the SAGB 
models using the range of values in Table 2 and baseline 
mixing parameters listed in Table 1. As an example of 
the rotation characteristics in the carbon core, Figure 11 
shows the angular momentum diffusion coefficients, spe¬ 
cific angular momentum, angular frequency, and mag¬ 
netic field profiles at first carbon ignition for an 8.0 M© 
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Model Number 

Fig. 12. — Angular momentum evolution of a 8 Mq star, with 
initial rotation (Q/£l cr it)i = 0.25 and overshoot f ov =0.016. The 
blue line shows the extent of the 4 He core, the black line the extent 
of the CO core and the black star marks the location of the first 
ignition. The white region at the top is where the star has lost 
mass. 
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Fig. 13.— HR diagram of two 8 Mq models: one non-rotating 
(blue line) and one with (Q/Q cr it)i ~ 0.2 (red) , for baseline mixing 
parameters. The diagram shows the evolution from H depletion to 
He depletion for both models. 

ZAMS model with (Q/Q crit )i= 0.2. Rotation is initial¬ 
ized by imposing a solid body rotation law at ZAMS. 
For (Q/fl crit )i= 0.2 this corresponds to an initial angular 
frequency of cu=1.522x 10 -5 rad s -1 and a total angular 
momentum of L=1.90xl0 51 erg s. At first carbon igni¬ 
tion, Figure 11 (middle panel) shows the carbon core has 
spun up a factor of ~ 30 to uo ~ 6xl0 -4 rad s -1 and 
rotates as a solid body. The largest angular momentum 
diffusion coefficient, by several orders of magnitude, at 
first carbon ignition is due to the Spruit-Tayler dynamo. 


The implied ~ 1 MG radial component of the magnetic 
field is shown in the bottom panel. DSI and SSI are not 
shown due to their negligible contributions in this model. 

Figure 12 illustrates the evolution of the specific an¬ 
gular momentum from ZAMS to first carbon ignition 
for the 8.0 M 0 ZAMS model initialized at ZAMS with 
(Q/Q cr it)i= 0.25 star and /ov=0 .016. During the main- 
sequence phase (model numbers 1000—1400) the specific 
angular momentum is uniformly distributed throughout 
the model. As the star ascends the RGB and then the 
AGB (model numbers 1400—3200) the specific angular 
momentum is extracted from the core and redistributed 
into the envelope, decreasing the specific angular momen¬ 
tum in the core by a factor of ~ 100. Boundaries between 
convective and non-convective regions are distinguished 
by sharp transitions in the specific angular momentum, 
with non-convective regions having the least specific an¬ 
gular momentum. The first ignition of carbon, off-center 
in this case, occurs around model 3700 and is marked. 

While Poelarends et al. (2008) claims that mass loss 
will strongly effect the final outcome of SAGB, we note 
that in our models (with a only one mass loss rate used) 
that mass loss seems to have a minimal effect on the star 
up to carbon ignition. In figure 12 the effect of mass 
loss is visible as the white region at the top of the fig¬ 
ure and only significant near the end of the core helium 
burning phase. While mass loss will extract angular mo¬ 
mentum from the star, the core’s angular momentum is 
unaffected. 

Rotation during the main sequence supplies a pro¬ 
longed source of hydrogen fuel that builds a slightly 
more massive helium core than non-rotating models (e.g., 
Maeder & Meynet 2000; Heger et al. 2000; Lagarde et al. 
2012). The increase in core mass (^ 0.05 — 0.2 M 0 ), 
shifts the effective temperature and luminosity. However, 
once core helium is depleted the two tracks converge on 
the AGB and there is little difference in the CO core 
masses between rotating and non-rotating stars. Figure 
13 shows a portion of the Hertzsprung-Russel diagram 
of two 8 Mq ZAMS models from H depletion to He de¬ 
pletion: one with (D/D cr ^)^=0.0 (blue curve), and one 
with {yt/VL cr i t )i= 0.2 (red curve). 

5.1. Evolution of Core Rotation 

Figure 14 (left panel) shows the evolution of the cen¬ 
tral rotation rate, Center, f° r different initial rotation 
rates for the 7.5 M 0 ZAMS mass model. During the 
main-sequence phase there is little evolution in the cen¬ 
tral rotation rate, spinning down by ~ 20%. Figure 14 
(left panel) shows the transition region between core hy¬ 
drogen burning to shell hydrogen burning and then to 
core helium burning. As a star leaves the core hydrogen 
burning phase (log 10 Age ^ 7.596 for the (D/D cr ^)i=0.1 
model), the nuclear energy in the core decreases. This 
causes a contraction of the core and subsequently a spin 
up of the core (Palacios et al. 2006). The convective core 
spins at a constant rate throughout the convective region, 
and spins faster than the outer-non convective layers. As 
the nuclear energy generated decreases, the convective 
region inside the core shrinks towards the center (Sills 
& Pinsonneault 2000). As it shrinks, magnetic fields in 
the outer radiative layers propagate inwards into regions 
which where previously convective. These fields act to 
slow down the core, removing the rotation differential 
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Fig. 14. — The evolution of the central rotation rate cn ce nter, for a 7.5 M@ at (Q/Q cr it)i=0.1,0.2,0.3,0.4,0.5. Left panel: Evolution from 
the TAMS to the ignition of core Helium burning. Right panel: Evolution during formation of the CO core and carbon ignition. Right 
panel insert: Zoom in on (r2/r2 cr ^fi=0.1,0.2,0.3 during the carbon ignition 


between the core and the envelope. 

Eventually the convective region recedes entirely al¬ 
lowing the magnetic fields to propagate through to the 
center, which causes the rapid spin down of uj center . At 
this point the core is still contracting, thus the core be¬ 
gins to spin up again. 

There are now two possible outcomes, seen in Figure 
14 (left panel), either the core has a second peak in 
core rotation ((fl/Q cr ^)i=0.1,0.3,0.5) or the core rota¬ 
tion plateaus ((fi/fi cr ^)i=0.2,0.4), the outcome of which 
depends on the sign of dcj/dM. Stars with dcj/dM < 0 
form a second peak in the core rotation profile. When 
convection restarts in the core, due to helium burning, 
the convective region expands outwards. As it does so it 
engulfs slower rotating material, which slows down the 
convective center of the star. Hence the second peak in 
core rotation occurs when convection restarts. For stars 
with dcj/dM > 0, a plateau is reached in the core ro¬ 
tation profile. Here, when the convection restarts and 
expands outwards it engulfs faster rotating material and 
thus spins up. However this forces the convective core 
to expand, thus as more material is engulfed in to the 
convective region, it will begin to slow the core down. 
Hence the maximum oj cente r occurs after convection has 
started. Typically, we find that the convective core will 
grow to ~ O.O5M 0 before the spin up ceases. While we 
believe these qualitative aspects hold, we caution against 
inferring quantitative predictions from this, as we have 
found that the sign of dcj/dM to be model resolution 
dependent. 

Figure 14 (right panel) show the evolution of cj cent er 
at the end of the core helium burning phase up to the 
start of carbon burning. Note the change in the scale of 
^center; the centers have spun down by « 25%, during 
the core helium burning phase. First, we see a glitch in 
the rotation rate at the end of the core helium phase. 
This is due to the same process that occurs at the end 
of core hydrogen burning. The convective core shrinks, 
allowing the magnetic fields to propagate inwards slowing 
the core down. This is countered by the core contracting 


and spinning up. We can also see that the initially faster 
rotating stars are evolving slower due to their ability to 
mix fresh fuel into the core (Maeder & Meynet 2000; 
Heger et al. 2000). 

We see that as the star forms its CO core the core spins 
up, due to the core contraction. Carbon ignition occurs 
at the glitches seen at log 10 cj C enter ~ —3.12 rad s -1 , 
when all stars have the same core rotation rate. These 
glitches occur due to the core expansion due to the car¬ 
bon burning events. Subsequent episodes of carbon burn¬ 
ing can be seen as smaller rotation glitches in figure 14 
(right panel insert). We can also see that the faster the 
star initially rotates, the later the ignition occurs, due to 
their slower evolution. 

5.2. Mass-Rotation Plane 

Figure 15 (top left) shows the location of the first 
carbon ignition in the mass-rotation plane for a fixed 
/ov— 0.016. Only models with 7 M 0 < M ZA ms < 8 M e 
feature off-center ignition, as models with Mzams < 

7 M 0 do not ignite carbon, while models with Mzams > 

8 M 0 feature central carbon ignition. Figure 15 (top 
right) shows the location where carbon flames and flashes 
are quenched. Only models having off-center ignition 
within the relatively narrow range 8.0 M 0 < Mza ms < 
8.2 M 0 does the flames or flashes reach the center. 

For a fixed ZAMS mass the ignition and quenching lo¬ 
cations are mostly independent of (0/O cr ^)^ values be¬ 
tween 0.0 and 0.5. This occurs because the transport of 
angular momentum transport from the core to the over- 
lying layers is efficient during the giant branch phases of 
evolution. Thus, regardless of the initial rotation rate, by 
the time the carbon core forms the central regions are ro¬ 
tating as a solid body with similar angular frequencies for 
a fixed ZAMS mass. Figure 15 (bottom left) shows the 
center angular frequency, Center, in the mass-rotation 
plane at the fixed baseline overshoot value. Note cj cent er 
only spans a factor of ~ 2 over the entire plane; all models 
rotate with similar angular frequencies at a fixed ZAMS 
mass. We can quantitatively explain, to first order, the 
rate at which the carbon core spins up between forma- 
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Fig. 15. — The ignition mass location (top left), minimum distance the carbon burning reaches to the core (top right), rotation of the 
center at ignition (bottom left) and CO core mass at ignition (bottom right) as a function of the initial ZAMS mass and initial rotation 
(Q/fi cr it)i at a fixed / ov = 0.016. White regions are models that do not ignite carbon. 


tion of the carbon core and first ignition of carbon from 
angular momentum conservation and the mass-radius re¬ 
lationship of polytropes. When the rotating carbon core 
forms, its total angular momentum is 

Li = IiUJi ~ aMiRfai , ( 20 ) 

where I is the moment of inertia, Mi is the mass of the 
carbon core, Ri is the radius of the carbon core, and q 
is a constant that depends on the density structure. At 
first ignition of carbon, the angular momentum of the 
more massive contracting CO core is 


L J = IfUf ~ CfMfRjujf . (21) 


Conserving angular momentum over this phase of evolu¬ 
tion gives 


ujf MiRj 

Vi ~ M f R) ' 


( 22 ) 


Assuming the non-rotating polytropic mass-radius rela¬ 


tion 

R ~ M (1 - n) / (3 - n) (23) 

applies at the first order, substitution into eq. 22 gives 

Uf _ Mf— 3w, )/ ( 3— ni) 

Ci ji M^~ ^ n /)/(3— n /) ? 

where rii is the polytropic index at the formation of the 
carbon core and n/ is the polytropic index at first car¬ 
bon ignition. For the angular frequency of the core to 
increase, ujf > uji , the polytropic index is restricted to 
be in the range 1 < n < 3. For example, for an 8 M 0 
ZAMS model with (Q/Ll cr i t )i= 0.2 at formation of the 
carbon core, we find Mi « 0.92 M 0 and rii ~ 1.5. Using 
a least squares fitting program to generate the polytropic 
index n for a sequence of MESA profiles between forma¬ 
tion of the carbon core and first carbon ignition (see Fig¬ 
ure 7 for an example), we find the left-hand side and the 
right-hand side of eq. 24 agree to within a factor of ^ 2 
for the 7 M 0 , 8 M 0 , and 9 M 0 models shown in Figure 
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Fig. 16. — The ignition mass (top left) as a function of the overshoot (/ov) and initial rotation at a fixed mass of 8 M©, minimum distance 
the flame reaches to the core (top right), rotation of of the center at ignition (bottom left) and finally the CO core mass at ignition (bottom 
right). White regions are models that do not ignite. Note the scale on the CO core mass is different to that of figures 15 & 17 


15. The center spins up, on average, by a factor of « 40 
between formation of the carbon core and first carbon 
ignition. Since the carbon core at first ignition rotates 
as a solid body with similar angular frequencies for a 
fixed ZAMS mass, the carbon core masses are nearly in 
the same state independent of the initial rotation rate, 
as shown in Figure 15 (bottom right). 

We find the carbon cores are rotating with periods be¬ 
tween 0.1 —1.0 days at first carbon ignition on the AGB. 
During the RGB phase of evolution we find the helium 
cores have periods of ~ 2.5 days, again independent of 
the initial rotation rate. Mosser et al. (2012) measured 
rotational splittings in a sample of about 300 red giants 
observed during more than two years with Kepler. They 
found these splittings are dominated by core rotation. 
Periods range between 10—100 days with larger periods 
for red clump stars compared to RGB stars. They in¬ 
ferred a ZAMS mass range of 1.2—1.5 M©, less massive 
than our rotating SAGB models. 

Stars with masses < 7M© will eventually form CO 
WDs after removing their outer envelopes. Between 


7 M© and < 8 M© where the carbon burning does not 
reach the core, our models suggest these stars will form 
CO+ONe hybrid WDs. Stars with masses > 8M© form 
ONe WDs as the carbon flames will burn away the 12 C. 
Electron capture supernovae are expected for stars with 
masses >9 M©, due to the CO core mass being greater 
than the Chandrasekhar mass (Eldridge & Tout 2004). 

5.3. Overshoot-Rotation Plane 

Figure 16 (top left) shows the location of the first car¬ 
bon ignition in the overshoot-rotation plane for a fixed 
ZAMS mass of 8 M©. For this case, overshoot is a domi¬ 
nant factor in setting the location of the first, off-center, 
carbon ignition. This first ignition of carbon can be made 
to occur at almost any mass coordinate within the car¬ 
bon core of the 8 M© model by varying the overshoot 
parameter. The f ov =Q .0 case, where convective mixing 
operates only within the Schwarzschild boundaries and 
does not extend beyond the MLT convective boundary, 
does not ignite carbon for any of the initial rotations 
rates. The smallest non-zero overshoot parameter in our 
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Fig. 17.— The ignition mass (top left) as a function of the initial ZAMS mass and overshoot (/ov) at a fixed (Q/Q C rit)i= 0.25, minimum 
distance the flame reaches to the core (top right), rotation of the center at ignition (bottom left) and finally the CO core mass at ignition 
(bottom right). White regions are models that do not ignite. 


grid, /ov=0.001, gives mass locations for the first igni¬ 
tion of carbon that are furthest from the center, closest 
to the outer boundary of the carbon core. Progressively 
larger values of the overshoot parameter generally move 
the location of the off-center ignition location closer to 
the center. 

Figure 16 (top right) shows the quenching location, 
where the flame and flashes die, in the overshoot-rotation 
plane for a fixed ZAMS mass of 8 M 0 . The flame and 
flashes approaches the center for nearly all the models; 
only models with / ov < 0.02 does the burning become 
quenched relatively far from the center (Denissenkov 
et al. 2013). Similar to our analysis the mass-rotation 
plane, Figure 16 (bottom left) shows carbon core rota¬ 
tion rate is approximately constant, to within a factor of 
« 2, regardless of the initial rotation rate. There is evi¬ 
dence for a weak dependence on the rotation rate to the 
overshoot parameter. Figure 16 (bottom right) shows 
the carbon core mass increases with increasing values of 
/ov, again nearly independent of the initial 
A larger core mass, in turn, leads to the first carbon 


ignition occurring deeper in the star. 

Comparing the results of the overshoot-rotation grid 
with the mass-rotation grid, we find the carbon core mass 
is the quantity that most strongly determines the struc¬ 
ture of the flame. For example, the boundary between 
cases that ignite off center and those that do not ignite 
carbon (ZAMS masses ~ 7 M 0 ) depends on whether the 
star can form a carbon core of ~ 1.05 M 0 , which is nec¬ 
essary to reach the critical density in equation 18. 

5.4. Mass-Overshoot Plane 

Figure 17 (top left) shows the mass location of carbon 
ignition in the mass-overshoot plane at a fixed ZAMS 
rotation of (fl/f} cr ^/=0.25. For / ov =0.0 the minimum 
mass needed to ignite carbon is 9 M 0 , « 2 M 0 greater 
than the baseline case, and models up to 11 M 0 ignite 
carbon off-center. While no overshoot may be unphysi¬ 
cal, even a small amount of overshoot moves the min¬ 
imum ZAMS mass for ignition considerably, down to 
7.8 M 0 . Increasing / ov decreases the required ZAMS 
mass to off-center carbon ignition and decreases the min- 
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imum ZAMS mass needed for central ignition of carbon. 

The width of the ZAMS mass range where the model 
stars ignite carbon off-center is approximately constant 
with respect to overshoot, AMzams/A/ ov ~ 1.6 M 0 . 
That is, overshoot uniformly moves the ZAMS mass 
boundaries where off-center carbon ignition occurs. For 
example, the sloped contours in Figure 17 (top left) show 
that when / ov =0.0, models in the mass range « 8.9 to 
11 M 0 have off-center ignition (as found by e.g. Siess 
2006, 2007) When / ov =0.008 this mass range shifts to « 
7.4 to 9.4 M 0 , and when / ov =0.016 the off-center carbon 
ignition range shifts to ~ 7.2 to 8.8 M 0 . 

Figure 17 (top right) shows the final fate of the car¬ 
bon burning flames and flashes, the quenching location 
in the mass-overshoot plane. With / ov =0.0, all flames 
and flashes reach the core (e.g., bottom right plot of Fig¬ 
ure 3). As the overshoot parameter increases, the carbon 
burning is less likely to reach the center. For / ov <0.01, 
the flame either has a single flash (similar to Figure 3 
top left panel for the 7 M 0 model) or undergoes a sin¬ 
gle flash then a steady state flame (similar to Figure 3 
middle right panel for the 8.2 M 0 model). Only when 
the overshoot parameter is large, / ov > 0.01, does an in¬ 
termediate evolution of a flash and then a steady state 
flame that does not reach the core (similar to Figure 3 
middle left panel for the 8.0 M 0 model). 

Figure 17 (bottom left) shows the central angular fre¬ 
quency at the point of first ignition in the mass-overshoot 
plane with a fixed ZAMS rotation of (Q/Q cr it)i= 0.25. 
We find that models which undergo central carbon ig¬ 
nition have lower angular frequencies than the off-center 
igniting stars. This is due to the carbon cores being spun 
up during the cooling phase, due to the contraction of the 
carbon core (see equation 24). As the heavier stars ignite 
carbon earlier they have less time in which to be spun up 
compared to the lower mass stars have at first ignition. 
As before, only models with carbon core masses greater 
than « 1.05 M 0 ignite carbon, as shown in in Figure 
17 (bottom right). As before, the maximum carbon core 
mass that ignites off-center carbon burning is ~ 1.2 M 0 , 
similar to that shown in Figures 15 and 16. 



Fig. 18.— ZAMS mass and final C/O core masses for non-igniting 
models in the mass-overshoot plane, measured once the 4 He shell 
has been depleted. 


Figure 18 shows the final CO core mass for the non¬ 
igniting models (white region in Figure 17). The maxi¬ 
mum mass for a CO core is 1.05 M 0 ; stars with heavier 
CO cores ignite carbon burning. We also find a trend for 
increasing overshoot to increase the final CO core mass, 
as noted previously. Doherty et al. (2015) reported a 
grid of models with the Monash stellar evolution code, 
over a range of metallicities, to investigate the fate of 
AGB and SAGB stars. Comparing our results with their 
solar metalicity results (their Fig. 6), we find that for 
a given ZAMS mass our rotating CO core masses are 
0.05—0.1 M 0 larger. Increasing the overshoot can mimic 
decreasing the metalicity in terms of the final CO WD 
mass. 



Fig. 19. — Temperature and nuclear energy generation rate pro¬ 
file of a carbon flame front within an 8 M 0 ZAMS model with 
(Q/£l cr it)i= 0.25 and f ov = 0.016. Red open circles show the mesh 
locations of the temperature while blue open circles show the ab¬ 
solute value of the nuclear energy generation rate. The distance 
between mesh locations within the body of the flame is <1 km, 
which is sufficient to accurately capture the nuclear burning and 
thermal transport processes. 


6. RESULTS FROM THE MIXING COEFFICIENT GRID 
STUDIES 

Table 5 shows overshoot has the most significant effect, 
on the location of the initial flame, with no ignition for no 
overshoot for the 8M 0 , (£l/£l cr it)i= 0.25 model, as long 
as the scale factor for the strength of angular momentum 
diffusion r] arn >0. As overshoot increases, the flame 
ignition occurs deeper in the star. This is due to changes 
in the 4 He and CO core masses during the stars evolution. 
For instance, comparing Figure 3 bottom left and bottom 
right, we can see the effect of overshoot for the 9 M 0 
model. Primarily, the model with overshoot ignites at the 
center, while the no overshoot model ignites off center. 
Without overshoot the CO core mass, the size of the 
helium shell and the ignition location are comparable to 
the 7 M 0 model with overshoot. 

There are two distinct populations in the thermoha- 
line models, those with small values of a t h, which ignite 
a flame at M « 0.5 M 0 and those with large a t h val¬ 
ues which ignite a flash at M ~ 0.8 M 0 (though only 
if / ov =0.001). Those models that ignite at M « 0.5 
M 0 , thermohaline mixing has little impact on the flame 
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TABLE 5 

Ignition locations in solar masses for the mixing grid. 


'Ham 


0.0 0.5 1.0 1.5 


Ot-sc 

OLth 

0.000 

Jov 

0.001 0.016 

0.020 

0.000 

Jov 

0.001 0.016 

0.020 

0.000 

0.001 

Jov 

0.016 

0.020 

0.000 

Jov 

0.001 0.016 

0.020 

0.000 

0.00 

a 


0.65 

0.28 


0.49 

0.15 

0.18 


0.43 

0.17 

0.16 


0.49 

0.17 

0.13 

0.000 

0.10 

0.50 

0.71 




0.46 

0.17 

0.05 

0.43 

0.42 

0.14 

0.05 


0.47 

0.18 

0.05 

0.000 

1.00 


0.63 

0.35 

0.38 

0.39 

0.46 

0.11 

0.18 


0.40 

0.16 

0.16 


0.48 

0.15 

0.16 

0.000 

10.0 

0.50 

0.65 


0.54 


0.49 

0.16 

0.15 


0.43 

0.18 

0.05 


0.43 

0.14 

0.16 

0.001 

0.00 

0.81 


0.65 



0.59 

0.21 

0.17 


0.83 

0.16 

0.14 


0.85 

0.23 

0.04 

0.001 

0.10 

0.79 

0.48 

0.28 

0.37 


0.83 

0.20 

0.05 


0.85 

0.14 

0.18 


0.69 

0.12 

0.05 

0.001 

1.00 

0.83 


0.58 

0.71 


0.84 

0.15 

0.17 


0.83 

0.18 

0.05 


0.71 

0.17 

0.15 

0.001 

10.0 

0.33 

0.69 

0.32 

0.49 


0.57 

0.18 

0.16 


0.62 

0.08 

0.05 


0.65 

0.11 

0.16 

0.010 

0.00 

0.19 

0.57 

0.45 

0.75 


0.63 

0.21 

0.16 


0.65 

0.17 

0.14 


0.63 

0.16 

0.05 

0.010 

0.10 



0.30 

0.39 


0.69 

0.20 

0.05 


0.62 

0.13 

0.04 


0.68 

0.12 

0.05 

0.010 

1.00 

0.73 



0.73 


0.66 

0.17 

0.18 


0.55 

0.20 b 

0.05 


0.85 

0.16 

0.14 

0.010 

10.0 

0.49 

0.57 

0.66 

0.36 


0.56 

0.18 

0.04 


0.63 

0.08 

0.17 


0.64 

0.10 

0.05 

0.100 

0.00 

0.28 


0.39 

0.76 


0.66 

0.12 

0.17 


0.84 

0.15 

0.14 


0.83 

0.16 

0.06 

0.100 

0.10 

0.54 


0.28 

0.37 


0.52 

0.20 

0.08 


0.58 

0.14 

0.17 


0.83 

0.11 

0.05 

0.100 

1.00 

0.51 


0.57 

0.73 


0.68 

0.15 

0.18 


0.83 

0.20 

0.06 


0.55 

0.17 

0.13 

0.100 

10.0 

0.69 


0.52 

0.55 


0.49 

0.18 

0.04 


0.58 

0.08 

0.15 


0.60 

0.12 

0.16 


a Ellipses represents models with no ignition 
b Baseline model 



Fig. 20.— Evolution of the timestep 5t (left y-axis) and Kelvin- 
Helmholtz thermal timescale r kh °f the carbon core (right y-axis) 
of an 8 Mq ZAMS model with (tl/Q cr it)i = 0.25 and f 0 v = 0.016. At 
model number « 2800 the timestep begins to decrease due to the 
increase in nuclear burning caused during core helium depletion. 
At carbon ignition, which occurs when the thermal timescale is 
in a local minimum at model number « 3700, the timestep is « 
10 years and decreases to w 1 year during the carbon flame and 
flashes. 

ignition location. However there is some variation due 
to thermohaline mixing before the flame ignites. Before 
the flame ignites, when e nuc >> there is a region in 
the vicinity of the ignition point that undergoes weak 
12 C+ 12 C burning (with e nuc << e v ). This weak burn¬ 
ing is able to drive a region of weak thermohaline mixing 
as a precursor to the vigorous carbon burning. Those 
that ignite at M ^ 0.8 ignite a flash, but under dif¬ 
ferent conditions to those that we predict for the other 
flashes. Here, there is thermohaline mixing between the 
CO and helium shell which allows us to form a small 
region (M « O.O5M 0 ) where X( 12 C) ~ X( 16 0). This 
higher fraction of carbon allows the ignition to occur at 
a lower density. This flash then prevents an ignition oc¬ 
curring deeper in the CO core where X( 12 C) « 0.3, which 


we assume in equation 14. 

Thermohaline mixing also effects the flame once burn¬ 
ing has commenced, the mixing pulls 12 C material from 
below the flame (Siess 2009). As the strength of thermo- 
haline mixing increases, the sub flashes, seen in figure 3, 
middle left panel, merge into one continuous flame, due 
to the increased carbon abundance. However as Siess 
(2009) showed this mixing eventually starves the flame 
of fuel preventing it from reaching the core. 

The effect of semiconvection is almost negligible, over 
the range of values considered here. Semiconvection acts 
near regions of convection, however it only acts for short 
periods of time in our models, thus has limited ability to 
change the composition of models before the formation 
of the CO core. In can however act during the carbon 
burning once the convective region has formed, mixing 
the burnt material with unburnt CO. Again this effect is 
small and plays a limited role the evolution of the flame. 
In figure 3, where we have secondary carbon flashes (top 
left, middle left and middle right panels), those flashes 
that occur near the 4 He shell can form brief semicon- 
vective regions across the shell and into the convective 
envelope. This may provide a way to detect the product 
of the flashes in the surface abundances, though these 
flashes (and the semi convective regions) are short lived, 
which will limit the material transferred. 

Changes in ? 7 am , which is a global scale factor on 
the strength of angular momentum mixing, primarily 
acts by changing the strength of thermohaline mixing. 
As r \ arn increases the amount of thermohaline mixing in¬ 
creases as well, which allows the mixing of material from 
the core into the flame (Siess 2009) to increase, though 
again this effect is small. 

Zaussinger & Spruit (2013) found for a 15M 0 model 
on the MS that the semiconvection mixing timescale is 
long (10 10 yrs) which explains why the semiconvection 
has little effect on these systems. Siess (2009) showed 
that thermohaline mixing has limited effect in the evo¬ 
lution up to the carbon ignition, due to the lack of the 
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Fig. 21. — Location of the ignition mass in the 7 Mq, 8 M©, and 9 M© models with a fixed / ov =0.016. Different values of spatial 
resolution (5 mes h(line style) and temporal resolution wt( line color) are show. Top left ( Q/Q cr it)i=0.0 , top right (Q/Q cr it)i= 0.25 and 
bottom (Q/Q cr it)i=0.5. Color/style combinations not shown do not ignite. 


3 He( 3 He,2p) 4 * & He reaction which is necessary to set up the 

mean molecular weight inversion needed for thermohaline 
mixing. Brown et al. (2013) propose a model for mixing 
by fingering convection in the parameter regime relevant 
for stellar (and planetary) interiors that is supported by 
three-dimensional direct numerical simulations. 

The angular momentum diffusion term has limited im¬ 
pact due to its implementation as a global scale factor 
on the angular momentum mixing process (Paxton et al. 
2013), thus the value itself is not physically motivated 
however we have varied to test whether missing sources of 
angular momentum, like internal gravity waves (Kumar 

& Quataert 1997), would have an impact. Given how¬ 
ever, that in Figure 11 we have shown that the individual 
diffusion coefficients can vary by 10 orders of magnitude, 
a change of « 50% in r\ arn is insignificant. For additional 
sources of angular momentum mixing to have an impact, 
they must be able to effect size of the CO core, like over¬ 
shoot does, to have a detectable impact. Compositional 
changes (thermohaline and semiconvection) are too weak 
to have appreciable impact on the ignition location due 


to their limited ability to change the CO core mass. 

7. RESULTS FROM THE SPATIAL AND TEMPORAL 
CONVERGENCE STUDIES 

Accurately capturing the nuclear burning and thermal 
transport within a steady-state, convectively bounded, 
carbon burning front, or within the time-dependent car¬ 
bon burning flashes, requires spatial resolutions < 2 km 
(Timmes et al. 1994; Ritossa et al. 1996; Garcfa-Berro 
et al. 1997). Siess (2006) use as many as « 50 grid points 
to describe the precursor flame between the convection 
region and the minimum in the luminosity profile be¬ 
low the convective region. Denissenkov et al. (2013) and 
Chen et al. (2014) use more than 100 mass zones. Fig¬ 
ure 19 shows the temperature and nuclear energy gen¬ 
eration rate profiles of a carbon burning flame in our 
8 Mq ZAMS model with (£t/£t C rit)i = 0.25 and / ov =0.016. 
Distances < 3700 km he ahead of the flame front, and 
distances between « 3700 km and ~ 3900 km contain 
the region where thermal conduction dominates nuclear 
burning. Distances between « 3900 km and ~ 4200 
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km contains the body of the flame front which reaches 
a peak temperature of « 7.5 xlO 8 K and peak energy 
generation rates of « 8.9xlO 6 erg g _1 s _1 . Distances 
> 4200 km contain the convectively bounded region of 
the flame. The critical temperature at which the heat¬ 
ing due to nuclear reactions equals the energy diffused 
away by neutrino and conductive processes in the steady 
state is about T cri t ^ 5.5xl0 8 K. The location of this 
critical temperature is marked in Figure 19 with a black 
cross. The profiles shown in Figure 19 also capture the 
flame structure with 1—2 km resolution with « 400 mesh 
points. The flame structure propagates inward toward 
the center at speeds of ~ 0.1 cm s _1 , consistent with the 
values reported in Timmes et al. (1994). 

Figure 20 shows the timestep and Kelvin-Helmholtz 
thermal timescale, r^h = GM%/R C L, for the CO core 
of the 8 M 0 ZAMS model with (Cl/Q crit )i= 0.25 and 
/ ov =0.016 model. As helium is depleted in the core at 
model number ~ 2800, the timestep begins to decrease 
from « 10,000 years due to the increase in nuclear burn¬ 
ing. At first carbon ignition, model « 3700, the timestep 
is « 10 years and decreases to w 1 year as the flame and 
flashes propagate towards the center. At model number 
~ 4500 the flashes have reached their closest approach to 
the center. The thermal timescale increases as the core 
increases in mass until the first ignition, where it then 
rapidly decreases due to the increased luminosity. The 
thermal timescale then peaks again shortly before the 
next ignition at model « 4400. This time however the 
flame generates less energy and the thermal timescale is 
reduced by a smaller amount, compared to the first igni¬ 
tion. On average the flame lifetime is ~ 10% that of the 
thermal timescale of the core. 

The location of first carbon ignition in the 7 M 0 , 8 M 0 , 
and 9 M 0 models as a function of spatial and temporal 
resolution is shown in Figure 21. Each model has the 
baseline f ov =0.016. Spatial resolution in MESA is gen¬ 
erally controlled by £ mes h, with smaller values providing 
an increase in the number of cells. Temporal resolution 
is loosely controlled by w t , the allowed change in the 
size of variables during a timestep, with smaller values 
decreasing the size of the timesteps taken. See (Paxton 
et al. 2011) for a detailed discussion of these two MESA 
control parameters. 

For zero rotation, Figure 21 (top left) shows all values 
of S mes h and w t give the same location of first carbon 
ignition, suggesting convergence has been attained. In¬ 
creasing the spatial resolution has little impact on the 
location of the flame, while increasing the temporal res¬ 
olution shows a slight decrease in the ignition location. 

For (£l/£lcrit)i= 0-25 (Figure 21 top right), the location 
of first carbon ignition depends on the values of S mes h and 
w t . For the 7M 0 case the highest resolution model (red- 
dashed), w t = 10 -5 and £ m esh=0.1, agree with our base¬ 
line model (green-dot-dashed), w t = 10 -4 and £ m esh=0.5, 
that there is no ignition occurs. At w t = 10 -3 , as the 
spatial resolution decreases the ignition point is pushed 
deeper into the star. At 8 M 0 , most of the models have 
converge around 0.1-0.2 M 0 ignition point, except the 
lowest resolution model with has a center ignition. Mod¬ 
els with w t — 10 _5 show little variation as £ mes h changes, 
while as w t increases in size the ^ mes hterm becomes more 
significant. These studies suggest our baseline values for 


^mesh and w t for off-center ignition are well within the 
convergence envelope. All models converge on a 9M 0 
star having center ignition. 

At (Q/Q C rit)i= 0-5 (Figure 21 bottom), for the 7 M 0 
case our baseline parameters agree with the highest res¬ 
olution model, for the lack of ignition. However, for all 
other values of w t and £ mes h there is considerable spread 
in ignition points. For the 8M 0 case the results have 
clustered around 0.05-2.0M 0 , except for w t = 10 -3 case, 
where the results have a spread of O.5M 0 . All models 
agree again the 9M 0 case has a central ignition. 

Overall, our baseline models agree within « O.1M 0 
of the highest resolution models we ran, for the igni¬ 
tion point. As the rotation rate of the star increases we 
can see that the choice of resolution terms becomes more 
significant and that there is a larger spread in possible 
values. Changing the temporal resolution has the most 
effect on the initial location of the flame. The choice 
of spatial resolution becomes more significant only as 
the temporal resolution decreases. Thus our choice of 
baseline parameters appear to be a good compromise in 
terms of precision of results and computational effort, de¬ 
creasing wt increases the computational time by a similar 
amount, while decreasing £ mes hincreases the memory re¬ 
quirements for the model. However, they also show a 
necessary requirement for carbon flame modelers to look 
critically at their choice of model resolution (Timmes 
et al. 1994; Ritossa et al. 1996; Siess 2006; Doherty et al. 
2010; Denissenkov et al. 2013; Chen et al. 2014). 

8. DISCUSSION 

We have investigated the detailed and global properties 
of carbon burning in SAGB stars with 2755 stellar evolu¬ 
tion models. These models consumed 200,000 core-hours 
(roughly 3 days per model) and yielded over 2 TB of dec¬ 
imated data (a limited number of MESA profiles were 
stored). To our knowledge this represents the largest 
block of compute resources used for a MESA survey to 
date. We note that every model ran from the pre-main- 
sequence to the end of carbon burning (if carbon ignition 
was achieved) without failure and without intervention. 

With these models, the location of first ignition 
whether off-center or central, the quenching location of 
the carbon burning flames and flashes, the angular fre¬ 
quency of the carbon core, and the carbon core mass have 
been surveyed as a function of the ZAMS mass, initial ro¬ 
tation rate, the magnitude of various mixing parameters 
such as convective overshoot, semiconvection, thermoha- 
line and angular momentum transport. We now compare 
our results to previous efforts and discuss methods for 
calibrating the f ov parameter within a given overshoot 
implementation. 

Georgy et al. (2013) found that rotation of a 9 M 0 
model can increase the lifetime spent on the MS com¬ 
pared to that of a non-rotating model. This increase in 
MS lifetime is caused by rotational mixing supplying a 
sustained amount of fresh hydrogen into the convective 
core. They include modifications to the stellar struc¬ 
ture equations due to centrifugal acceleration described 
by Kippenhahn et al. (1970); Endal & Sofia (1976), as¬ 
suming the angular velocity is constant on isobars (Zahn 
1992). Georgy et al. (2013) also adopt an instantaneous 
method of overshoot with d over /H p = 0.10 applied to 
the H- and He- burning boundaries. For a non-rotating 
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9 M 0 model at the end of core He burning they find a 
convective core to total mass ratio of M cc /M tot « 0.10 
while our corresponding rotating model yields a larger 
value of Mcc/Mtot ~ 0.15. We find a more modest differ¬ 
ence between the non-rotating and rotating 8 M 0 model 
of Mcc/Mtot ~ 0.14 and M cc /M tot & 0.15, respectively, 
as shown in Figure 13. 

While our MESA models use a similar implementa¬ 
tion for rotation, our calculations differ from Georgy 
et al. (2013) in that we include the effects magnetic 
torques which aid in significantly inhibiting the spin up 
of the convective core of the star during its evolution. 
For example, we find for an 8 M 0 ZAMS model with 
(£l/£l C rit)i ~ 0-5, an angular velocity at the center of 
the core of logio ^center ~ —3.4 rad s -1 at the start 
of carbon ignition, compared to a more rapid value of 
logio ^center ~ —0.9 rad s _1 when internal magnetic 
fields are neglected. Magnetic field torques that inhibit 
spin up of the stellar model results in less massive con¬ 
vective cores due to the less efficient rotational mixing. 
We also find that magnetic torques can account for the 
less drastic shift in luminosity on the HRD for an 8.0 M 0 
ZAMS model with (Q/Q crit )i^ 0.2 (See Figure 13), con¬ 
trary to the larger differences in the rotating model HRD 
tracks shown in Georgy et al. (2012). 

For their stars past 2DU, Doherty et al. (2015) found 
the CO core mass for a 7 M 0 ZAMS model to be 
^0.8 M 0 and increases to 1.375 M 0 for their 9.5 M 0 
ZAMS models. In mild contrast, our models predict the 
CO core mass for a 7 M 0 ZAMS model to be 1.05 M 0 
and the highest mass star to produce a Chandrasekhar 
core to be a 9 M 0 ZAMS, assuming an overshoot of 
f ov =0.016. These differences are likely due to the treat¬ 
ment of the convective boundaries, with Doherty et al. 
(2015) using a search for convective neutrality rather 
than a convective-decay prescription, leading to differ¬ 
ences in the size of the 4 He and CO core masses. 

Siess (2006) found for models without overshoot that 
the ZAMS mass range which ignite carbon off-center 
is 9—11.3 M 0 . This is comparable to our / ov = 0.0, 
rotating models (Figure 17), which yield a value of 
M up « 8.8 M 0 and M mas > 11 M 0 . 

Denissenkov et al. (2013) found for a 9.5 M 0 ZAMS 
mass with no overshoot or thermohaline mixing, an off- 
center ignition mass of 0.665 M 0 with the flame pro¬ 
ceeding to the center, consistent with our results. With 
/ov =0.007 they found carbon ignites off-center but the 
flames and flashes do not reach the center. In contrast, 
we find in this case the model star undergoes a central ig¬ 
nition. We speculate this difference is due to Denissenkov 
et al. (2013) only including overshoot once the CO core 
has formed, where we include it from the pre-MS on¬ 
wards. Thus, the CO core in Denissenkov et al. (2013), 
/ov=0.007 model will be smaller than our models, and 
hence the ignition occurring off-center. 

While our models are not completely comparable to 
those of Jones et al. (2013), who use / ov =0.014 except at 
the base of a burning region where they use / ov =0.007, 
they find a 8.2 M 0 model ignites off-center while models 
with M zams > 8.8 M 0 centrally ignite. We find that for 
models with 0.007 < f ov < 0.014 that a 8.2 M 0 will ignite 
off-center while only models with M zams > 9.4 M 0 will 
always centrally ignite. 


Arguably the biggest uncertainty in stellar models is 
the treatment of convection. The overshooting param¬ 
eter in particular, regardless of how it is implemented 
within a specific numerical instrument, critically influ¬ 
ences all outputs of stellar evolution (e.g. Maeder 1975, 
1976). Figure 17 in particular demonstrates that the 
properties of carbon burning in SAGB models is not an 
exception, especially the range of ZAMS masses that ex¬ 
perience off-center ignition. Testing on a small number 
of models suggests that the most significant location for 
overshoot is in regions of He burning, followed by car¬ 
bon burning. Regions with H burning or no burning 
show little difference in ignition location with respect to 
changes in overshoot. The effect of convective overshoot 
on the stellar models considered in this work are in agree¬ 
ment with previous work by Siess (2007) who showed 
that for / ov = 0.016 applied at the edge of the convec¬ 
tive boundary, M up can transition from 8.90T 0.10 M 0 to 
7.25=L 0.25 M 0 for Z = Z@. We find a similar transition 
where M up « 8.8M 0 for our rotating, / ov = 0.0 model, 
which shifts to a value of M up ~ 7.2 M 0 for / ov = 0.016 
(See Figure 17). Gil-Pons et al. (2007) found similar re¬ 
sults upon investigating a grid of zero metallicity stars 
with / ov = 0.12 using an instantaneous overshooting 
formalism (Herwig et al. 1997), contrary to the diffu¬ 
sive approach used in this work. They find a value of 
M up « 6.0 M 0 and M mas « 7.8 M 0 . The adoption of 
instantaneous overshooting, as well as Z Zq are likely 
to contribute to the modest discrepancy in values of M up 
and M mas . 

Traditionally the value of the overshooting parameter 
for a given overshooting model is calibrated by fitting 
isochrones against the width of the terminal age main 
sequence in color-magnitude diagrams, or the surface 
abundances, of young and intermediate age clusters (e.g., 
Maeder 1976; Maeder & Mermilliod 1981; Mermilliod & 
Maeder 1986; Schaller et al. 1992; Herwig 2000; Vanden- 
Berg et al. 2006; Kamath et al. 2012). Photometry and 
spectroscopy of binary systems offer another avenue for 
calibration of overshooting because these measurements 
can provide the radii, effective temperatures, and masses. 
In addition both components of the binary need to lie 
on the same isochrone and fit their respective evolution¬ 
ary tracks (Schroder et al. 1997; Pols et al. 1997; Ribas 
et al. 2000; Claret 2007; Meng & Zhang 2014; Stancliffe 
et al. 2015) High-precision high-cadence space photom¬ 
etry from the CoRoT and Kepler missions opens up a 
newer method for calibration of overshooting and other 
mixing processes in stellar interiors (Neiner et al. 2012; 
Montalban et al. 2013; Tkachenko et al. 2014; Guenther 
et al. 2014; Aerts 2015). 

MESA implements the time-dependent treatment of 
convective overshoot mixing of Herwig (2000) with the 
traditional calibration method leading to / ov =0.016. It 
is unknown if this value of / ov in this specific overshoot 
model applies to masses other than the ones used for 
calibration, is consistent with values derived from binary 
systems or asteroseismology, or if it applies to advanced 
burning stages of stellar evolution. However, we have 
shown that for a dense grid of SAGB models taken to the 
end of carbon burning, utilizing our adopted baseline pa¬ 
rameters, values of M up rsj 7.0M 0 and M mas - 8.4 M 0 
are nearly independent of initial of ZAMS rotational val¬ 
ues of (D/f l cr it)i ~ 0.0 - 0.5. While our SAGB models 
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have been evolved from the pre-MS phase through the 
end of carbon burning, for models whom do not ignite 
carbon, and those that ignite carbon off-center, the ini¬ 
tial rotational rate may play a larger role in the final 
rotational rates of the WD that will eventually be born. 
For a given ZAMS mass and overshoot parameterization, 
we suggest that strong claims of carbon burning quench¬ 
ing at an appreciable distance from the center to yield 
hybrid CO + ONeNa white dwarfs should be viewed with 
caution. 
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